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I. INTRODUCTION 


Quantum Electrodynamics (QED) is an extremely successful theory and its predictions 
have been conhrmed experimentally with outstanding precision. This success has corrobo¬ 
rated the general approach originally developed by Feynman, Schwinger, and Tomonaga, and 
based on the renormalization theory (see, e.g., the monographs Feynman, Schwinger, 


and Tomonaga’s approach was extended by Furry p to the description of processes occur¬ 
ring in the presence of some classes of classical, background electromagnetic helds (these 
classes exclude, for example, those electromagnetic helds, which render the QED vacuum 
unstable under electron-positron pair production, see [6|). The application of the so-called 
Furry picture is required when the interaction of electrons (and positrons) with the external 
held has to be taken into account exactly. The corresponding mathematical and physical 
conditions depend on the structure of the background held. Here, we are interested in pro¬ 
cesses occurring in the presence of a strong laser beam and the controlling parameter is 
the so-called classical nonlinearity parameter ^ = \e\Fo/mu}o, where Fq is the amplitude of 
the laser wave and cjq is its central angular frequency (m and e < 0 are the electron mass 
and charge, respectively, and units with c = 1 are employed throughout) |7l4lfll|. If ^ > 1, 
it is necessary to include the ehects of the laser held exactly in the calculations from the 
beginning. Since a realistic laser held has a complex space-time structure, this is in gen¬ 
eral an impossible task as it requires the exact analytical solution of the Dirac equation in 


3y 

m- 


the presence of the background laser held. However, the task has been accomplishec 
Volkov in the ideal case of a laser held approximated by a plane wave [ll[ (see also 
Starting from the seminal paper 12|, the corresponding electron wave functions (Volkov 
states) have been extensively employed to investigate various strong-held QED processes in 
a plane-wave held from the hrst-order ones as nonlinear Compton scattering 13l-l27l| , nonlin¬ 


ear Breit-Wheeler pair production 


13 


15 


28 


35| . nonlinear Bethe-Heitler pair production 


36l-l42|. and electron-positron annihilation [^, to higher-order processes like nonlinear Dou¬ 
ble Compton scattering 44-46] and trident pair production j^, 48|. In order to calculate a 
second-order process like nonlinear double Compton scattering, the Green’s function of the 
Dirac equation in the presence of a plane wave (Volkov propagator) has to be used [7|, [sl. 


49- 


5l[ |. The Volkov propagator is also a useful tool to inves 


presence of a laser held like the polarization operator 


igate radiative corrections in the 


52| (see also 53j). It is remarkable. 
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however, that originally the polarization operator has been determined independently in 


by employing the operator technique, which avoids t 
The same operator technique has been employed in 


le direct use of the Volkov propagator. 


55| to calculate the mass operator in a 


plane wave and in 5^ to investigate the process of photon splitting in a plane wave. The 
representation of the mass operator and of the polarization operator found in 55| and in 


54| via the operator technique has been later widely exploited to study radiative corrections 


and vacuum-polarization effects in plane-wave helds 


57 


Breit-Wheeler pair production via the optical theorem 


6^2| and the total rate of nonlinear 


63 | (see also 64|). 


We emphasize that nowadays the parameter region ^ ^ 1, where the laser held has to be 
taken into account exactly in determining the electron dynamics, can be routinely entered 
experimentally by means of optical laser beams. In fact, the threshold .^ = 1 at an optical 
laser photon energy of 1 eV requires a laser intensity of the order of 10^® W/cm^ and the 
present record intensity achieved in the laboratory exceeds this value by about four orders 
of magnitude |65|. On the other hand, we point out that the condition ^ > 1 is not related 
to the importance of pure quantum effects like the photon recoil in the emission of radiation 
or the creation of an electron-positron pair in laser-photon collision. The importance of such 
effects is determined by the ratio of the laser held amplitude Fq and the the so-called critical 


held of QED: Fcr = nF/h\e\ = 1.3 x 10^® V/cm = 4.4 x 10^^ G { 2 , 6,l6^. The ratio Fq/F^ 


m 


the laboratory frame is not a Lorentz invariant quantity and therefore it cannot represent by 
itself the physical parameter controlling the importance of quantum ehects. In the case of 
processes primed by an electron with initial four-momentum and colliding with a plane- 
wave held, it turns out that the important parameter is the so-called quantum nonlinearity 
parameter y = {{kop)/muo){Fo/Fcr) where hk^ is the laser photons four-momentum. 

The parameter y represents the amplitude of the laser held in units of the critical held of 
QED in the initial rest-frame of the electron. Analogously, in the case of nonlinear Breit- 
Wheeler pair production by a photon with four-momentum q^, the quantum nonlinearity 


parameter is k = {{koq)/muJo){Fo/Fcr) [7-ll0|. which corresponds to the amplitude of the 


plane-wave held in units of Fcr in the hnal center-of-momentum frame of the created electron 
and positron. Entering the strong-held QED regime where not only ^ > 1 but also y ^ 1 
(here, we limit to processes primed by an electron), is a challenging task nowadays as the 
amplitude of the critical held of QED largely exceeds available held strengths (the critical 
held amplitude corresponds to a laser intensity of the order of 10^® W/cm^). Thus, even 
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considering npcoming 10-PW facilities as APOLLON-lOP 67|] (expected intensities of the 
order of 10^^ W/cm^) or future lasers like the Extreme Light Infrastructure (ELI) ^ and 
the Exawatt Center for Extreme Light Studies (XCELS) 69(] (expected intensities of the 
order of W/cm^), this regime can be entered only by employing ultra-relativistic 

electrons. In this respect, it is also worth pointing out that the expected intensities are 
assumed to be achieved by spatially focusing the laser energy almost down to the diffraction 
limit. Unavoidably, also temporal compression of the laser energy into few-cycle pulses is 
required but Volkov states and propagators can still be applied for temporally compressed 
plane waves. On the contrary, they cannot be employed if the background held is a spatially 
focused laser beam. We mention here that the so-far unique experimental campaign on 
strong-laser held QED [to, 71| has been carried out with a laser of wavelength of the order 
of 1 pm and a relatively large pulse spot area (~ 60 pm^), such that the experimental results 
were well reproduced theoretically within the plane-wave held approximation. 

The requirement of tightly focusing the laser energy also in space in order to reach exper¬ 
imentally ultra-high intensities has motivated us to investigate the possibility of determining 
analytically the electron wave functions for tightly focused laser helds under certain realistic 
approximations (^we have mentioned, the problem in general does not admit an exact ana¬ 
lytical solution) 721]. In 72||, we have assumed that the electron is ultra-relativistic and that 
its initial energy e is the largest dynamical energy in the problem 73||. We emphasize that 


the formalism employed in 


72| and here below does not apply only to the electron but rather 


to any charged particle, provided that the dynamical conditions of validity of the approach 
are fulhlled (implicitly one also requires, of course, that the particle can be assumed to be 
stable for the whole duration of the interaction). In this respect, although below we will 
refer for dehniteness to the electron and the wave functions will contain the electron charge 
and mass, it should be borne in mind that the results can also be employed in principle for 
other charged particles. Now, the meaning of the assumption on the initial energy of the 
electron being the largest dynamical energy in the problem was to ensure that the electron 
is barely deviated from its initial direction, which was chosen as the reference direction 
for introducing light-cone coordinates. In the experimentally relevant case of an electron 
initially counterpropagating with respect to a tightly focused laser beam, the assumption 
implied in particular that although the laser amplitude is such that ^ 1, as it is envisaged 

in planned experiment on strong-held QED, the transverse momentum of the electron in 
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the laser field (of the order of m^) still remains much smaller than e. Under this approx¬ 
imation, we constructed the electron wave functions within the Wentzel-Kramers-Brillouin 
(WKB) approximation, by including terms scaling as 1/e j^, 75|. We notice that the above 
assumption e S> well fits present and near-future experimental conditions envisaged to 
test strong-field QED with intense lasers. For example, even at the next-generation facility 
APOLLON-lOP the quantity will be of the order of 80 MeV 0 and, in order to reach 
the threshold x = 1 where nonlinear QED effects become essential, electron energies of 
about 500 MeV are required. Moreover, electron beams with energies of about 4 GeV have 
been recently demonstrated experimentally with laser-plasma accelerators 76| (conventional 
accelerators provide electron beams with energies even exceeding 200 GeV, see, e.g., [77|). It 
is worth pointing out that the results in [72| and those presented here are also valid if ^ ~ 1. 
Now, at ~ 1 the formation region of strong-field QED processes as nonlinear Gompton 
scattering and nonlinear Breit-Wheeler pair production in a plane-wave field is of the order 
of the laser wavelength js] and the tight focusing of the laser field is expected to essentially 
alter the probabilities with respect to the plane-wave results. On the contrary, in the regime 
^ S> 1 the formation region of the mentioned strong-field QED processes in a plane wave has 
been found to be ^ times smaller than the laser wavelength |8|. In this regime, one would 
expect, for example, that the total probability of the mentioned strong-field QED processes 
in a tightly focused laser pulse can be obtained by averaging the probability calculated in 
the presence of a constant crossed field, after replacing the latter with the local value of 
the tightly focused laser field (local const ant-crossed field approximation) js]. However, it 
has been found recently that the local const ant-crossed field approximation in some circum¬ 
stances does not reproduce accurately the harmonic structures of radiation emission spectra 
in nonlinear Gompton scattering 78|]. Although the differences might be explained as a 
consequence of interference effects among different emission points such interference 

effects are non-local and this also calls for a thorough analytical investigation of strong-field 
QED processes in the presence of a tightly focused laser beam. Moreover, such an anal¬ 
ysis can also shed light on the accuracy of the local constant-crossed field approximation 
at intermediate values of the parameter ^ and for those fields where the local value of the 
parameter ^ varies from large values to values of the order of or less than unity. 

In the present paper, we continue and extend the investigation started in 72| in the fol¬ 
lowing directions. First, in Section |TT] we present a simpler expression of the wave functions 
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found in 


72| . which is equivalent under the approximations employed. Moreover, we gener¬ 


alize the previous results allowing for a longitudinal direction of light-cone coordinates non 
necessarily coinciding with the initial momentum of the electron. Second, in Section IIIIl we 
construct the corresponding scalar and spinor propagator under the same approximations. 
We note that a general expression of the WKB propagator for a scalar particle has been 


obtained in 


80|, by applying the so-called path-dependent formulation of a gauge theory 


based on the path-integral representation of the propagator. Here, we employ the operator 


technique, which we have already applied, e.g., in 


81 


to calculate the electron quasiclas- 


sical wave functions and propagator in the presence of a combined atomic and plane-wave 
held. After determining the propagator, we show that the found wave functions can be 
obtained by applying a limiting procedure to the propagator itself. Finally, in Section HVl 
the main conclusions of the paper and possible future investigations are presented. 


II. ULTRA-RELATIVISTIC WKB ELECTRON WAVE FUNCTIONS IN A 

FOCUSED LASER FIELD 


We consider an ultrarelativistic electron moving in a background electromagnetic held, 
described by the four-vector potential A^{x) in the Lorentz gauge = 0. The four-vector 

= {e,p), with e = \/m? + indicates the asymptotic four-momentum of the electron, 
i.e., either in the remote past or in the remote future. Even though the considerations 
below apply to electromagnetic helds of rather general space-time structure, we have in 
mind the case where the background electromagnetic held represents an intense, few-cycle, 
and tightly focused laser beam. Thus, we also assume that the held tensor = 

d^A^{x) — d'^A>^{x) = {E{x), B{x)) is localized on a space and a time region of a few 
cubic typical laser wavelengths and a few typical laser periods, respectively, and that it 
has a maximum amplitude Fq. The typical laser wavelength and period correspond to a 
typical laser angular frequency, which is indicated as cuq and which is such that the classical 
nonlinearity parameter ^ satishes the strong inequalities: m <C £ (see 72|]). As we 

have already hinted in the Introduction, under the above assumptions the electron is only 
slightly dehected from its asymptotic direction by the background held in the physically 
relevant situation where it is initially counterpropagating with respect to the laser held and, 
from this point of view, the regime ~ 1 can also be investigated with the present method. 
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Thus, it is convenient to introduce the light-cone coordinates 
t + n ■ X 


T = 


x± = X — [n ■ x)n, 


= t 


n ■ X 


( 1 ) 


for a space-time point x with space-time coordinates = {t, x) and to indicate the four 
light-cone coordinates via the corresponding capital letter X\ = (T, a;x,0). In [^, for 
the sake of simplicity, we have chosen n to coincide with p/\p\ as this essentially simplihes 
the solution of the classical equations of motion. However, as it will be clear below, this 
assumption is not strictly necessary and more general expressions of the electron wave func- 


< 

r\j 


<C P+, where 


tions can be obtained by assuming that the direction n is such that |p+| 
p+ = p — {n ■ p)n and p+ = [e + n ■ p)/2 are three of the four light-cone coordinates of 
the asymptotic four-momentum. The on-shell condition implies that the remain¬ 

ing light-cone coordinate p_ = e — n ■ p can be expressed in terms of the other three as 
P- = (m^ -|-p+)/2p+. The light-cone coordinates of an arbitrary four-vector = (n°, v) are 

dehned in the same way as for the four-momentum: n+ = +n-v)/2, v_[_ = v — {n-v)n, and 

n_ = -n-v. The same definition is extended to the Dirac gamma matrices 7 ^ = ( 7 °, 7 ), 

which are intended in the Dirac representation: 7 + = ( 7 ° -|- n ■ 7 )/ 2 , 7+ = 7 — (n • 7)n, and 
7 _ = 7 ° — n ■ 7. Moreover, the derivatives 


A 

&T 


d 


d 


A 

d(j) 




d 


d{n ■ x] 


( 2 ) 


( 3 ) 


dt d{ri ■ xy 

with respect to the light-cone coordinates T and 0 can be derived from the relations 

d 1 d d d 1 d d 

dt 2dT~^ d(j)' d{n ■ x) 2dT 90’ 

whereas the derivatives with respect to the transverse coordinates form the two-dimensional 

vector V+ = d/dx±. Also, the scalar product between two four-vectors and can be 

written as 

{uv) = M+n_ -h M_n+ — u_i_ ■ v_i_ (4) 

and the four-divergence of a vector held V^{x) as 


91A dV 
^ dT ^ d6 


V, -y, 


( 5 ) 


Below, for the sake of completeness, we summarize the important steps we followed in 


72| in order to solve the Dirac equation 


[Yiihd^ - eAp) - m]i/> = 0, 


( 6 ) 









for specific asymptotic conditions on the electron wave function 'il){x\A) and within the 
parameter regime specihed above. Recall that in terms of the light-cone coordinates the 
largest dynamical energy scale of the problem is , which is approximately equal to e in the 
ultra-relativistic regime. Now, based on the general argument that the de Broglie wavelength 


of an ultra-relativistic particle is very small, we applied the WKB method 75(] anc 
or a solution of the Dirac equation of the form '0(a:|y4) = exp[iS'(x|R)/h]<p(a:|y4) 


we looked 


Sal-y. In 


72l |. we have observed that, under the assumption that m <C e, the WKB method can 


be applied if the condition h/wQS 1 is fulhlled, with wq being the laser waist size. This 
condition is safely satished in the case of a tightly focused, optical {wq ~ 1 pm) laser held, as 
h/e ~ mXc/s Ac, with Ac = h/m = 3.9 x 10“^^ cm being the Compton wavelength. By 
writing the electron wave function as 'iIj{x\A) = exp[i5'(a:|y4)/h]<p(x|R), the Dirac equation 
becomes 


+ eA^) -h m]ip = 


( 7 ) 


By hrst neglecting the term proportional to h, one obtains 


[y^{df,S + eAfj) +m](p = 0 


( 8 ) 


for the bispinor ^p{x\A) at zero-order in h. This equation admits a non-vanishing solution 
only if det[7'^(9^S' -|- eA^) -|- m] = 0, which implies that the quantity S'(a;|74) has to satisfy 
the Hamilton-Jacobi equation 


{d^S + eA^,){d^^S + eA^^) - = 0 


( 9 ) 


and that it can be identihed with the classical action 


Q. 


In 


72|, starting from the solution 


of the classical equations of motion, we employed the method of characteristics in order to 
construct the action S{x\A). Here, we follow a different procedure, which allows us to easily 
obtain a more general expression of the action, including the possibility that p± does not 
vanish as it was in [^. First, we pass to light-cone coordinates and we initially assume that 
is the electron four-momentum at a hxed reference “time” Tq, which will be set to ±oo 
later. Thus, we write the action as S'(XIH) = Sp{X;Tq\A) = —{p+cj) + P-T — p_i_ ■ a;_i_) -|- 
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6Sp{X;Tq\A). The quantity 6Sp{X;Tq\A) satisfies the equation 


P+ 


+ eA_ ) + p_L ■ {'V±5Sp - eAj 


d6Sp 


, ddSp , dSSp 


dSSpd6Sp 1 2 2a A 1 2 42 


± 


eA± ■ 'V±6Sp + p- 


d6Sp 

d(j) 


+ eA_|_ I — 0. 


( 10 ) 


This equation can now be solved iteratively by setting 5S'p(X; TqIA) = (5S'p°^(X; TqIA) + 
5S'p^^(X;To|y4) + • • •, with (5S'p"'^(X;To|y4) = 0{l/p\). It is clear that at the leading order 
we can write 


5S^^\X-To\A) = -e / dT'/l_(T', 0), 


( 11 ) 


'To 


where the same initial conditions as in [7^ have been assumed for the four-vector potential, 
i.e., y4^(To, a;_i_, 0) = 0. By substituting this expression of 5S'p°^(X; To|y4) into Eq. flTUD . it is 
easy to show that the first-order correction TqIA) reads 

1 r'^ 

)i'?r-T„i4'i = — / 


6S^^\X-,To\A) = — dr 
P+ Jto 


ep± ■ A±{r, x±, r To) - -e‘^A\{T', x±, 0; Tq) 


( 12 ) 


with 


A^ix-,r) = A^{x) + v^ / drA_ir,x^,(j)). 


(13) 


'To 


In order to relate the results below to the finding in 72|], we observe that the quantity 
Aj_(X; To) can be written in terms of the electromagnetic field as 

rT 


A^{X;r) = - / dr[E^{r,x^,cl)) + nxB^ir,x^,cl>)], 

JTo 


(14) 


which is equal to —Gp(r, ro) in the notation employed in |72|. Also, we notice that up to 


first order in l/p+, the dependence on (j) can be neglected in the external four-potential. 
In fact, as it can be seen in Eq. (3) in |72|, it is 0(T) — 0o = 0(l/p+)- Correspondingly, 
Eq. (6) there shows that, by neglecting the derivatives of the four-potential with respect 
to (j), the conjugated momentum p_|_ — eA+(T, a;_i_) is conserved up to terms 0(l/p+) with 
respect to the leading one. Since here we will neglect the dependence on 0 of the four-vector 
potential (one can assume for simplicity that 0o = 0), it is convenient to introduce the 
three-dimensional light-cone coordinates X = (T, x±) and define the effective transverse 
vector potential 

A^{X-To) = A^{X) + V^[ dT'A_{T\x^). (15) 

JTo 
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Thus, we can write the action Sp{X-, Tq\A) up to hrst order terms in l/p+ as (for the sake of 
notational simplicity and since no confusion can arise, from now on we do not write explicitly 
the upper index (1)) 

rT 


SpiX;To\A)=-{p+cP + p.T-p^-x^)-e / dT'A^T^x^) 


'To 


P+ Jto 


dT' 


epx ■ -4±(r', a:±; To) - -e^^i(r', T„) 


( 16 ) 


which generalizes the expression found in 72|. It is also worth observing that the present 
approach allows to easily hnd higher-order corrections in l/p+ in the action by means of Eq. 

OB. 

After determining the action, the construction of the first-order electron wave function 
ippp iXlA) = To|A), with a being the spin quantum number at Tq, proceeds as in 

72l | and the final result is 

^/jpAX;To\A) =e^sAX;TolA)/n(^^ e . A±ir,xr,To) 


V 2 Jrp^ P + 

tS ■ [B{T\ x±)-nx E{T', 


u 


'p,a 




= ‘le 


(17) 

B. 


where S = —i 7 ^ 7 ^ 7 ^ 7 , Up^a is the usual constant free bispinor normalized as u'p^^Up^cr 
and a unit quantization volume has been assumed. Apart from the more general expression 
of the action Sp{X;Tq\A), this expression of the wave function is equivalent to the one found 


m 


72|. 


Now, as it is observed, e.g., in j^, when performing the matrix elements of basic strong- 
held QED processes like nonlinear Compton scattering and nonlinear Breit-Wheeler pair pro¬ 
duction in the ultra-relativistic regime, the contribution of those terms in the pre-exponent 
proportional to the unit matrix is suppressed with respect to those having a different matrix 
structure. The reason is that, since all particles involved in such processes almost propagate 
along the same direction and thus the photon polarization vectors e^, with A = 1, 2, almost 
lie on the plane perpendicular to this direction, the matrix element of the interaction matrix 
7 ■ ex between two constant bispinors vanishes at leading order in l/p+. As a result, the 
leading-order terms proportional to the unit matrix and the terms 0(l/p+) with a different 
matrix structure in the pre-exponent in Eq. flTTll give contributions of the same order of 
magnitude to the transition matrix elements. Thus, in order to calculate probabilities of 
strong-held QED processes at the leading order within our approximations, all hrst-order 
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corrections in the pre-exponent proportional to the unit matrix can be neglected. This is 


also formally corroborated by observing that, from the derivation presented in 


721, the pre¬ 


exponential function in Eq. flTTll proportional to the unit matrix results from the overall 
exponential factor exp (l/2p+) dT'd^d^Sp{T', a;j_, cj); To|y4) (see Eq. (11) in the Supple¬ 
mental Material in 72|), which is proportional to the well-known Van Vleck determinant 861 


(see also 


871. l88| l and which can be neglected within our approximations, as correspond¬ 


ing to a quantum correction to the action. The above argumentation allows us to write 
the wave function 'ijjp^f^{X]TQ\A) in a simpler form than in Eq. (IT7|) but equivalent within 
our approximations. In fact, by exploiting the approximated relation j^Up^a ~ n ■ jUp^a in 
the hrst-order term in the pre-exponent in Eq. flTTl) . it can be easily shown that the wave 
function 'ijjp^(j{X;To\A) can also be written as 


'4^pAX;To\A) = Ep{X;To\A) 


u. 


'p,cr 


^/2i’ 


(18) 


where the matrix 


Ep{X;To\A) = e^SAx-,To\A)/h 


1 - 


2p+ 


7+7_l • A±{X;To 


(19) 


is the analogous of the Ritus matrix in the case of a focused laser held (see B)- 

This form of the wave function is very suggestive. In fact, apart from the inconsequential 


phase factor exp 


— i{e/h) dT'A_(T', x^) (see also Eq. flT6|) ). which does not contribute 
when evaluating transition probabilities of strong-held QED processes, the structure of the 
state 'ipp^aiX] Tq\A) is the same as that of a Volkov wave function for a plane wave propagating 
along the direction —n , and described by the four-vector potential A'^{T) = (0, Ax(T),0), 
with the substitution 

A^{T) ^ A^{X-To). (20) 


This turns out to be especially interesting, by calculating the loca 
turn nonlinearity parameter x starting from its general dehnition 


value x{X) of the quan- 


x{X} = ^yiEE(X) + p X B(X)f - Ip . E(X)Y. (21) 

In fact, up to the leading order in l/p+, one obtains that x(V) is independent of 0 and it 
is given by 


xiX) 


p+ 1 

dA^{X-To) 

P+ 1 

dAAX) ^ , 

m Ecr 

dT 

m Ecr 



( 22 ) 
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which exactly corresponds to the local value of the parameter y for a “plane wave” with 
four-vector potential A^{X;Tq) = (0, A_l(X;Tq), 0) (note that the laser phase corresponds 
to 2T). This is in agreement with the general observation that an ultra-relativistic charged 
particle “sees” an arbitrary electromagnetic field locally as a plane wave [^. Of course, the 
same correspondence also occurs by introducing a local definition E(X) of the parameter ^ 
as the ratio between x(^) and the quantity h{kp)/rn? (see, e.g., [5^) or by starting from 


the manifestly Lorentz and gauge invariant definition of the parameter ^ proposed in 


18|- 


A related observation is in order here. In fact, in general, a Dirac bispinor 'ip{x\A), which is 
solution of the Dirac equation (E)), transforms as 'iIj{x\A) '0'(a:|A') = eyi]){—ieg{x))'il>{x\A) 
under the gauge transformation A^^{x) A'^{x) = A^{x) + d^g{x). Now, let us consider the 

gauge transformation, which makes the component A'_{X) of the transformed four-vector 
potential vanish, such that the latter resembles the four-vector potential in the plane-wave 
case (recall that the component Aj^{X) does not play any role under our approximations). 
The corresponding gauge function can be chosen as 5 f(X; Tq) = — dT'A_{T', x^) and the 

transformed wave function resulting from the one in Eq. flTSl) would exactly correspond to 
a Volkov wave function in the plane-wave potential A±{T) with the substitution A±{T) —)■ 
AY(X;To)(=^x(X;To)). 


For the sake of completeness, we conclude this Section by writing the in- and out-states 
with positive and negative energy, which can be obtained from the general expressions in 
Eqs. fll6l) and flT8|) by means of the limiting procedures To —)■ ±cxo, as mentioned in 72 |. 


By employing the same notation as in [2[ for the negative-energy constant bispinors M_p,_cr 


normalized as u^_p _^u-p-a = 2e, the final result is: 




At.(X\A) 

i,^°;%{x\A) 


^ ^iS^A^{X\A)/h 
^ ^iS^:f{X\A)/h 


2p+ 

e 


2p+ 


7+7±- 

7+7i-A“’(V 


u 


'p,a 


v^’ 


U 


'p,a 




1 + TT—7+7± ■ 


2p+ 


XL—r\ — 


p-G 


1 + 


2p+ 


7+71 ■ xr\x 


ATe ' 

ti—p^ — Q- 

JTIT’ 


(23) 

(24) 
(26) 
(26) 





















13 


where 


^g^(X|A) = T(p+0 + P-T -p±x^)-e I dT'A.{T\ x^) 


P+ J-c 


dT' 


epx • A^r\T', x^) T x^) 


S’^;^\X\A) = t{p+4> + P-T ~P±-x^) + e dT'A_{T', x^) 


P+ Jt 


dT' 


epx . '(T'.xx) T ^e^Ar''(T',a:^) 


2 * (out), 2 , rpl 


Arm = A^{X) + Vj. / dT'A_{T\ *±), 

J —oo 

/ oo 

dT’A.{T',x_L)^ 


( 27 ) 


(28) 

(29) 

(30) 


III. ULTRA-RELATIVISTIC WKB PROPAGATOR IN A FOCUSED LASER 

FIELD 


In the present section we derive the electron propagator nnder 


Section mi and we closely follow the approach developed in 


81 


he same assnmptions of 
, based on the operator 


techniqne. The electron propagator G{x,x'\A) satisfies the eqnation 

riihd^ — eAm)) ~ x'\A) = 5{x — x'). 

It is convenient to write the electron propagator in the form 


(31) 


G{x, x'\A) = — eAm)) + x'\A) (32) 

snch that the “sqnare” propagator D{x,x'\A) satisfies the eqnation 

{ihd^ — ~ eA^{x)) — {x) — 

where = (1/2) [ 7 ^, Before carrying out the spinor case, one takes advantage in con¬ 
sidering the simpler scalar case. By indicating as (x, x'\A) the corresponding propagator, 
one has to solve the equation 

— eA^((r))(ihcl^ — 624 ^( 0 :)) — D^'^\x,x'\A) = 5{x — x'). (34) 

In order to work within the operator-technique approach, we introduce the momenta oper¬ 
ators corresponding to the light-cone coordinates as 

= -ihV±, (35) 


D(x, x'\A) = 6{x — x'), (33) 
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and we express the propagator via the Schwinger proper-time representation 

• poo 

D^°\X,X'\A) = -- / (lse^W^'>^‘2P^^T-Pl+2e{A+PT+A_P^+A^-P^)WA^-m^]^(^j^ _ 

^ Jo 

Here, the four-dimensional (5-function 6{X — X') can be expressed as 

and the four-vector can be identified with the four-momentum flowing through the prop¬ 
agator, which also explains the use of the symbol in accordance to the results in Section 
im By replacing this expression of the (5-function into the propagator, we obtain 


(5(X -X') = J 


D^°\X,X'\A) = - 


h 


ds 


d^p . 


{27rh)' 


g ilp+{fl>-<l}')+P-{T-T')-p±-{xj_-x'^_)]/h^i{p'^-m^)s/h 


X exp \ i- [2P^Pt - PI- 2(p+Pt + + P± ■ Pi 


+ 2e[H+(Pr — p_) -|- A-{P^ — p+) -|- A± ■ (P± + P±)] + 

(38) 

Analogously to what we have already carried out in the Section [TTl we will apply the prop¬ 
agator to the investigation of processes, where the particles are ultra-relativistic and whose 
propagation direction is only barely altered by the laser held. Thus, we will keep only terms 
in the exponent up to order l/p+ and we will neglect all quantum corrections. Since, as we 
will see below, the relevant proper-times s scale as l/p+, we can already neglect all terms 
containing the quantity p_ in the last exponential in Eq. fl38|) . We point out that, even 
though the propagator contains an integration over all possible four-momenta p^, the con¬ 
servation laws and the dynamics in the processes of interest here, will automatically single 
out the relevant integration regions in the four-momentum and in those regions it will be 
P- |p_l| <C p+. Analogously to the procedure followed in [Sll, l82|, this also allows us to 
neglect the held component A AX) in Eq. (|38i) . which, after neglecting the terms propor¬ 


tional to p_, always appears as a correction to the large quantity p_|_ in the combination 
— eA^{X) (we will come back later on this point). Finally, as we have already pointed 
out in Section m we can also neglect the dependence of the background four-vector poten¬ 
tial on the light-cone variable 0. These considerations allow to substantially simplify the 
propagator D'P\X^X'\A)^ which can be written in the form 


D^^\X,X'\A) = - 


h 


ds 


^ P ^-i\p+(<l>-4>')+P-(T-T')-p±-(xj_-x\)]/h^i(p^-7n'^)s/h 


i2Tihy 


(39) 


X e 


-i{s/h){Pl+2{p+PT+p±-P±)+2e[A-P+-A^-{P^+pjyi+e^A\} 
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This expression of the propagator precisely corresponds to the starting expression of the 


propagator in 


81 


Now, the possibility of neglecting the quantity p_ in the last ex¬ 


ponential brings an additional simplification as the integration in p_ is trivial in Eq. fIMll . 
Indeed, it essentially provides a h-function (5(T —T' —2sp_|_), which in turn allows to perform 
the proper time integration: 


D^^\X,X'\A) = 


k 


^ f d^P. 


27ih 




(27rh) 


'-L )/h^-i(p'^, +m'^)s/h 


(40) 


X e 


-i{so/h){Pl+2{p+PT+p±-P±)+2e[A-P+-Aj_-{Pj_+pj_)]+e‘^A‘]_} 


where Sq = {T — T')/2p_|_. It could have been formally more convenient to perform the 
integration p_ later and to avoid the appearance of the T-dependent quantity Sq in the 
exponential, which does not commute with the operator Pt- On the one hand, this does not 
introduce ambiguities under our approximations, as all quantum corrections in the exponen¬ 
tials in Eq. fHU]) are neglected. On the other hand, we have ensured that the final results 
are unchanged if one had followed the alternative procedure and perform the integration in 
p_ later. Thus, we can ignore this subtlety here below. 

Now, the task is to disentangle the operator 




(41) 


where s is a numerical parameter. First, it is convenient to disentangle the operator 
exp[—2i{s/h){ p+Pt + p± ■ -Pl)], which corresponds to the free motion of the particle. Thus, 
we write 

0{s) = (42) 

with 

— ^-i{s/h){Pl+2{p+PT+p±-P±)+2e[A-p+-Aj_-{Pj_+pj_)]+e^A\}^2i{s/h){p+PT+p±-P±) 

The operator L{s) satishes the differential equation 

{Pl + M^-is)p+ - ^±(s) ■ iP±+P±)] +e^A\{s)]] , (44) 

where we have introduced the short notation 


y4_(s) = A_{T — 2sp+, x±_ — 2sj?_l), 
Ax(s) = A±(T - 2sp+,xx - 2spx). 


(45) 

(46) 
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Since we keep only leading-order terms in h, we can neglect the non-commutativity of the 
operators at different values of s and write 




(47) 


Here, we have exploited the fact that the operator exp[—2i(s/h)(p+PT + P± ■ -Pl)] in 0{s) 
acts on the unit function once the operator 0{s) itself is substituted back into X'\A). 

Now, since s scales as l/p+, the only operator in the exponent of 0{s) of the order of u nity 
is the one containing the quantity p+ ds'A_(s'). According to the procedure in 
it is convenient to disentangle all other operators by writing 

0(s) = 


81 


(48) 


where 


_g-(*/^‘)/o <is'{^l+2elA_(s')p+-Aj^(s')-(rj^+pj^)]+e^Al(s')]j 
X ghAl/o* ds'[Pl-2eA^{s')-{P^+pi_)+e^A\{s')\ ^ 


(49) 


By applying the same procedure as in the case of the operator L(s), we can write the operator 
M(s) under the same approximations as 


and 


_g-(2*P+/?i)/g® ds'A-{T-2s'p+,x^-2s'p^+2e ds" Aj^{s'')-2s'P±^ 

0(s) =e“(2®P+A)/o’ ds'A-(^T-2s'p+,x±-2s'p±+2efg ds" Aj^{s")-2s'Pj^ 
X g-hAi/o* ds'[Pl-2eA^{s')-{P^+p^^+e^A\{s')\ 


(50) 


(51) 


Now, since the argument of the second exponential already scales as l/p+, we can neglect the 
non-commutativity of the components of the operator P± with the corresponding coordinates 
and approximate the operator 0(s) as 


0(s) ds' A-{T- 2 s'p+,x^- 2 s'p^+ 2 e ds” A^{s'')- 2 s'P^^ 

X ds'[ 2 eA^{s')-p^-e^A\{s')\ 

rygOAl/o" ds'[2eA^{s'yp^-e^A\{s')\ 

X g-(2*P+A)/o'* (r-2s'p+,a;x-2s'px+2e/Q“ ds”Aj^(s'')- 2 s'Pj^ 


(52) 


We observe that the operator Pl cannot be neglected in the argument of the held component 
A-(X) because the corresponding operator in the exponential is of the order of unity. By 
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replacing this expression of the operator 0{s) in Eq. (|40D and by performing the shift 
p± —!■ p± + (x \ —x'l )/2sn (recall that from classical electrodynamics the quantity \x_[_ — x'^\ 
scales as l/p+, see (72|), we obtain 

i f dp+ _ /ft f (Ppi 


2h J 2Tih |p+| 




4.0 


j {2'RhY 

X !q° ds{2eA^{s)-[p^+{x^-x'^)/2sQ\-e^A\(s)} 

^ g-(2*P+/?i)/o*° dsA-iT-2sp+,x^--l^{x^-x'^)-2sp^+2e ds' A^[s')-2sP^^ 


(53) 


where from now on the meaning of the abbreviation Ax(s) is 

' ' ' -2sp_L 


A_l(s) = AAT-2sp+,x±_ - {x^-x'^ 

So 


Aj_ (T—2sp+, x_[_ - ix_}^—x\ 

So 


(54) 


The last approximation is justihed because all terms containing Ax(s) in Eq. fl5^ scale 
already as l/p+ and because the integration in p \ will produce either quantum corrections 
(due to the factor exp(—isoP^/^)) IS2|) or corrections scaling as l/p+ (due to 

the factor exp [/(2e/h)p_L • (isAj_(s)]). 

At this point, the next step is to shift the variable p± as p± p± + (e/so) /q*° dsA±(s). 
Since after this shift the variable p± will only produce quantum corrections, it can be ne¬ 
glected in the argument of the held component A_{X) and the remaining Gaussian integral 
in p± can be performed analytically. The result is 


D^°\X,X'\A) = - 


Airti? / 2'Kh 


exp — —m So 
|T-T'| h 


X exp 


hsn 


X± — X, 


PSQ 1 2 

+ e (isA_i_(s) 

Jo 


— [ dsA']_{s) 


h 


X e 


-(2*p+/rt)/o*° dsA- (T-2sp+,x^-j^{x^-x'^)+2eJ^ ds'A^{s')-2e-^ dsA±{s)-2sP±^ 


(55) 

As we have already mentioned, our aim here is the determination of the propagator within the 
WKB approximation. Thus, we have neglected the transverse quantum huctuations brought 


about by the original presence of the term proportional to p± in A_(X) 


81 


and, apart 


rom this aspect, the expression in Eq. fl55|) is in agreement with the corresponding one in 
8ll. l82| . By expanding the last exponential operator up to terms of the order of l/p+ and 


by updating the abbreviation A_(s) according to Eq. fl5Tl) . i.e., 

A_(s) =A_(t 


2sp+,x_i - {x_i 

So 


X 


(56) 



























we obtain 
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D^^\X,X’\A) = 


AtiK^ J 27ih 
i 


exp 


|T-T'| 


h 


m So 


X exp 


hs 


0 


x±_ — x\_ 


rso 


n 2 


+ e / dsA±(s) 


j fSO 

a2 

h 


dsA, ('s') 


>0 


X exp 


4ie^p+ 


'•so 


(is[9j_y4_(s)] 


r-so 


I ds'A±(s') -/ dsA±(s) 

L-'o 'So Jo 


X e 


— {2iep^lh) dsA_(s)+(4iep+//t) Jq^ ds s[5xj4-(s)]-Pj^ 


(57) 


where we have introduced the two-dimensional gradient d± with respect to the whole trans¬ 
verse vector variable of the held to which it is applied (the transverse vector variable 
Xj_ — {s/so){x± — x'j_) in Eq. (I57|) h The hnal task is to disentangle the last exponential in 
Eq. fl57|) . It is convenient to introduce a parameter p and to write 


R[p) = _ j^^^p\^^{Aiep+/h) fP dss[d^A-{s)]-P^ ^ 


with 


N{p) = dsA-{s)+{iiep+/h) dss[S^^_(s)]-Pxg-(4iep+/fi)/p'" ds s[ax^-(s)]'-PL 


(59) 


By differentiating with respect to the variable p and by noticing that in general [P±, A_(s)] = 
—ih{l — s/so)d±A_{s) (see Eq. fl56|) h we obtain that 


N{p) = e“(2iep+/;i)/o^ds|A_(s)+4ep+(l-^) 


[axA-(s)]./;ds'y[axA_(s')]} 


(60) 


In conclusion, since the operator on the right of N{p) in Eq. fl58|) acts on the unit function 
(once the operator R{p) is substituted back into Eq. ([57])), we obtain 


D^°\X,X'\A) = 


dp+ - 


ip+{<f>-tl)')/h_ 


AtvIt? / 27rh 


6 ^( 50 ) 
\T -T'\ 


exp 


* 2 

-m So 
n 


2iep+ 

h 


cso 


dsA_{s) 


X exp 

X exp 

X exp 


z 

hso 


rsQ 


-f e / dsA±{s) — / (isA^(s) 


r'so 


'0 


'■so 


Aie^p^ 

Jo 

%ie^p\ 


(is[^_L74_(s)] • 


[ ds'A±{s') - [ dsA±{s) 

Jo So Jo 


h 


ds ( 1 —^ ) [9 _l74_(s)] ■ I ds' s'[^_l74_(s')] 


( 61 ) 
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As a check on this expression of the propagator, we go back to Eq. (IM|) and we notice that 
the propagator D^^\x,x'\A) equivalently solves the equation 


[{ihd'ij^ + eAn{x')){ihd'’^ + eA^{x)) — D^^\x,x'\A) = 6{x — x'). (62) 


This can be easily proven by writing the propagator D^^\x,x'\A) in the operator form as 
D^^\x,x'\A) = {x\D^^'>{A)\x') and by noticing that Eqs. fl3il) and fl6^ correspond to the 
operator equations 


[{P,-eA,(X)){P>‘-eA‘‘{X))-m^]D^'>\A) 

(63) 

= D^^\A) [{P^ - eA,{X)){P>^ - eA>^{X)) - m^] = 1, 


where X^ and P^ are four-position and the four-momentum operator, respectively (with 
an abuse of notation we have employed the symbol X^ only in this equation to indicate 
the four-position operator, whereas in the rest of the paper it indicates the generic space- 
time point in light-cone coordinates). Now, Eq. fl6^ clearly shows that the propagator 
D^^\x,x'\A) has to fulhll the general symmetry relation D^^\x,x'\A) = D^^\x',x\ — A). In 
order to verify that indeed the propagator in Eq. fl6T|) satishes this symmetry relation, it 
is convenient to perform the changes of variable: a) —p+, such that the quantity Sq 

remains unchanged (after the exchange X -H- X'), b) s —)■ — s -|- Sq (s' —)■ —s' -|- Sq) in the 
integrations in s (s'). 


A more explicitly symmetric form of the propagator D^^\X, X'\A) can be obtained by 
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performing the change of variable t = T — 2sp+ {t' = T — 2s'p^) in the integrals in s (s'): 


D^^\X,X'\A) = - 


Airh? \T -V 


exp 


——e / dTA_{T) + —e 


i x± — X 


h 


IT' 


h T-T' 


dTAAr) 


IT' 


X 


Q-^p+ 

27rh 


T-T' 

2p+ 


exp 


^ T -T' ip+ {xt - 
h 2p+ h 2 T — T' 


X exp 


X exp 


i e 


h2p+ IT-T' 

,2 rT 


r A 1 

2 .T ^ 

/ drA_L(r) 

- / drAlA) 

_Jt' 

Jt' j 


% e 


hr[^_L^-(r)] 


T' 

T-t 

T-T' 


r T' - T 

i - JTTf 


dr'A At' 


X exp 


/ i TT' 
\h2p+T -T' 


X 


IT' 


fT' 


dr (l - [^±2l_(r)] ■ dr' [ax2l_(r')] 

dr (^1 - dr' (^1 - [^x^_(r')]|y 


where 


2l_(r) = A_{t,xt- - x'^j , , 

T-t 

Aj_(t) = Ai ( t. Xx - - ai'i 


(64) 

(65) 

( 66 ) 


Now, let us consider the exponential function 

T-t 


£:(X, X') = exp y (ir[a_L44_(r)]- 


T-T 


dr'A At') - 


'T' 


T -T 
T-T 


dr'A At') 


(67) 


in Eq. (16^ . where we notationally reminded that a dependence on the transverse variables 
Xt and x'^ is brought about by the components of the four-vector potential (see Eqs. fl65|) - 
fl66|) h By integrating by parts, it can easily be shown that the quantity £{X,X') can be 
written as 


S{X,X') =exp 


i e 
hp+ 


2 rT 


dr A 1 (t) 


iT' 


AAx(Ti X,X')- AtT L drAATr, X, X') 


( 68 ) 
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where 


AAx(r;X,X') 


T 


T _T> 


T — t' 

dT'd±A_ ( r', x±_ - —— 


T' 


/ dr'd^A- ( r', a;_L — 

T-T' Jt' V T - T' 


T -T' 
T-r' 


(ail - 


(69) 


By manipulating in a analogous way the last exponential function in Eq. fIMll . one can 
hnally write the propagator X'|a4) in the compact form 


D^^\X,X'\A) = 


1 


47rh2 |T - T' 


exp 


i / , . / X z a;_L — I 

.-ej^^drA_{r) + -e^^ 


dr A I ('r'l 


IT' 


X 


„ /T - T' 


27rh 


g-ip+{<t>-<t>')/HQ 


X exp 


/ 


\ h2p+ I T-T' 


2p+ 

fT 


exp 


——m 


T — T' f (a;_i_ — a; 


/ 


+ 


_ h 2p+ h 2 T-T 


f dTAT{T;X,X') 

Jt' 

- £dTAl{T;X,X’)^ 


Here, we have introduced the effective transverse vector potential 


(70) 


ATr, X,X’) = A^{t) + Ax4x(t; X,X'), 


(71) 


which exactly corresponds to the analogous effective transverse vector potential At{X]Tq) 


introduced in Section [ITl in the sense that it reduces to and to {t; x'j_ 

in the limits T' —)■ —cx) and T ^ ex, respectively (see also the discussion below Eq. flTHli ). 
In this respect, we want to show that, analogously to what we have observed in Section 
m the propagator can be written in an alternative, equivalent form such that, after per¬ 
forming the gauge transformation A^{x) —?■ A'^{x) = A^{x) -|- d^g{x) with gauge function 
g{X]To) = — dTA_{T, x±), the resulting propagator in the new gauge coincides with the 
scalar Volkov propagator with the plane-wave transverse vector potential A±(T) being re¬ 
placed by A'^(X; To) = A±{X; Tq) (see, e.g., [49|, ISOj). Although this alternative expression 
of the propagator is formally asymmetric with respect to, e.g., the transverse variables x± 
and x'^, it is instructive to sketch the proof of the above statement by reminding that under 
our approximations the quantity \xt — x'^ \ scales as l/p+. By starting from this observation. 


1 (out) / 
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it is easy to show that 
— e f dTA_{T) + e 


JT’ 

—e 


x± — x\ 
T -T' 


dr A I(r) 


/T' 

r-T' 


dTA_{T, x^ + e 


'To 


j dTA_(T, x'^)+e ^^ _ ■ j dTAjT,xj_;Ti 


(72) 


Moreover, since the argument of the held-dependent exponential inside the integral in in 
Eq. (1701) is already proportional to l/p+, one can easily see from Eq. fl6^ that 


A±iT-,X,X') ^A±{T,xr,To) 


T -T' 


pTo rj^l pTo 

J dTV±A_{T', ^±) - J. J dTV^A_{T',x^). 


(73) 


T -T 

The proof of the above statement is complete once one observes that, since the difference 
A^{t-X,X')-A^{t, ^uTq) does not depend on r, then 


T-T 


-\ r dTA^{T-,x,x')\ - r dTAi{T-,x,x' 

Ut' j Jt' 


T-T 


rT 1 2 

dTA±{T, x±;To) 


UT' 


(74) 


drA^ir, x±]To 


IT' 


Coming back to Eq. fITOD . we make two additional remarks. First, the argument of the 
exponential function outside the integral in coincides within our approximations with 
the quantity —i{e/h){AXfj_/AT) dTA'^{T), with = (AT, Ax±, A0) = (T — T', x± — 
x'j_,(j) — (j)'), as A0 scales as l/p+. Second, the quantity p+ appears in the held-dependent 
part of the integrand in Eq. fITOl) only in the denominator, which can be used as an a 
posteriori qualitative indication that the held component A+(X) can be neglected in the 
propagator within our approximations (see the discussion below Eq. fl38|) h 

Now, the integral in is Eq. (TfOll can be performed analytically. We hrst remove the 
0-function by splitting the integral from —oo to -|-oo into two integrals, one from —oo to 0 
and the other from 0 to -|-cxo. Then, by performing the change of variable = yAT/2 in 
both integrals and by combining them, we obtain 

1 


D^^\X,X'\A) = - 




exp 




. i Ax 
dTA_{T) + —e- 


X 


dyexp ( —— 


T' 

i fAX^ 


h AT 


dr A j (t) 


IT' 


h\ 4 


y+ \m^ + 6m^(X,X')] - 

y 


(75) 


where the quantity 

W(Jf.X') = 


f 1 r 9 

r 1 r ,1 

'1 

[ATiA^^TX.X-)- 

[at L 



(76) 
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is the local generaliza 


focused laser field 


49 


ion of the mass correction due to the laser held to the case of a 
l^see also 


90 


9l|h Finally, by performing the change of variable 


1 / —)■ l/l/, it is clear that the integrand coincides with that of the free (and of the Volkov) 
propagator, with the substitution —?■ M^{X, X') = m? + 5m?{X^ X') > such that 
the propagator can be written in a very compact form as 92(] 


D^^\X,X'\A) = - 


exp 


X 


AttK^ 


i AX, 




h AT 


IT' 


dTA^{T) 


1 M(X,X') 




17] 


:Hf*(Af(X,X')VAA'2 -iij) 


(77) 


( 0\ 

where H) (z) is the Hankel function [93| and ?7 is a positive inhnitesimal quantity. 

Now, by introducing the reduced propagator D^^^X, X';p+\A) such that (see Eq. flTOll ) 


D(°)(X,X'|A) = / ^e-^P+^^-'^''^/^D^^\X,X';p+\A), 
I 27ih 


(78) 


it can easily be shown that it generates the in- and out-states with positive and negative 
energies = exp [iS'£p'^°“*^(X|A)/h] = exp(=Fip+0)$|(.p'^°'^*^(X|A) of a scalar 

particle in agreement with the electron wave functions found in Section 2 and according to 
exactly the limiting procedure and substitution rules as in Eqs. (68)-(73) in [8^ (one has 
only to remember that here the in- and out-states are dehned as those states which reduce 
to the free states in the remote past and future, respectively, also for the negative-energy 
states): 


lim Cl“>(X,X';p+|7l) = 

T-foo 
Itc^l—>-00 


lim D^^\X,X']p+\A) = 

T'^-oo 
|cc'j^|—>-oo 

lim D^^\X,X']-\p+\\A) = 

T'^oo 
|aj^ 1^00 

hm D'"l(V,X';-|p+||7l) = 

T —>--00 

|£Cj^ |—>00 


1 


47rh2 T P 

1 pi{P±-'m'^)\T'\/2hp+ 


Ann? |T'| 

/2h\p+\ 


4i?>(X|/l) 


1 


47rh^ 


T' 




(out) 
-q ( 


X|A) 


1 A{Q±-^'^)\T\/^fi\P+\ . 

—«f5*(X'|A) 


AttK^ 


where 


p'" =P+{ 1, 1 




m 


2 T2 2p2 


X 




X 


22-/5 


m 


|T| 

= \p+\ 

q'^ = b+| 


1 — 1 - 
’ 'j'/ ’ 

1 1 

, |T|’ 


X 


/2 


V+ > 0, 
P+ > 0, 

P+ < 0, 

P+ < 0, 

2 

m 


2T'- 


X-, 


2pi 


m 


2P 2p\ 


(79) 

(80) 
(81) 
(82) 

(83) 

(84) 
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Now, by looking at Eq. (13^ for the square propagator D{X^X'\A)^ we see that the 
correction in the exponential operators induced by the new spin terms scales as l/p+ and 
it is independent of h (which means that it is a first-order quantum correction in h, as the 
exponent in the propagator always contains the overall factor 1/h). Thus, by introducing 
the corresponding reduced propagator D{X, X';p^\A) such that 


D{X,X^\A) = 
it is easy to show that 


X'; p+\A), 

27rh 


(85) 


DiX,X';p+\A) = |l + • [Ax(X;To 


Ax(X';To)]U(o)(X,X';p+|kl), (86) 


where it should be noticed that the difference A.±{X;Tq) — A.±{X';Tq) does not actually 
depend on Tq. Also, we point out that, analogously to the procedure followed in Section UTl 
in deriving the matrix structure in Eq. fl86|) we have retained only terms up to l/p+, we 
have neglected corrections proportional to l/p+ to the unit matrix, and we have exploited 
the fact that at the leading order in l/p+ the second term in the curly brackets in Eq. fl86ll 
is always multiplied by a matrix of the form 'y ■ ex, with ex being the polarization vector of 
a photon (almost) propagating along n (see Section [IT]), and then by a free bispinor. 

Now, the spinor propagator G{X, X'\A) in light-cone coordinates can be obtained from 
Eq. (l3^ . By introducing the reduced electron propagator G{X, X';p+|A) such that 


G{X,X'\A)= f ^e-*P+('^-'^')/''G(X,X';p+|A), 
./ znn 


(87) 


also in the spinor case, it is 

G(X,X';p+|A) = 


7"^ +7 

h 7 _l ■ + eA_i_{X)) + m 


( 88 ) 


D{X,X'-p+\A). 


By exploiting the arguments below Eq. fITT]) and by noticing that 
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ifiA _ eA.(X] 
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X exp 


—i 
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'To 
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= exp 


—z 
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dTA_{T,x_i_) + i— / dTA_{T,x'j_) 
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( 89 ) 


7^*^^+7 p'^+ 7±-{ihV± + eA±{X;To))+m 
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it can be seen that also in the spinor case the propagator can be written alternatively in an 
equivalent form such that, after performing the gauge transformation mentioned below Eq. 
CB. the resulting propagator in the new gauge coincides with the spinor Volkov propagator 
with the plane-wave transverse vector potential Aj_(T) being replaced by A'j_(X;To) = 
.4._l(X;To). In order to show this statement, it is also important to observe that within 
our approximations and after the commutation in Eq. (1891) has been carried out, the action 
of the operator in Eq. flHHj) on the background four-vector potential can be ignored. 

This formal equivalence indicates that the propagator G(V, X'\A) can also be written in the 
compact form 


G(X,X'\A) = 


E,(X-,T„\A) , . E,(X'-,T„\A), 


(27rh)"^ — m'^ + irj 


(90) 


like in the plane-wave case js], where p = where in general Ep{X] TqI^) = To|A) 7 °, 

and where the dependence on Tq hnally compensates out. 

Finally, it can be easily shown that the analogous limiting procedures as in Eqs. (74)-(77) 
in j^l allow to obtain the in- and out-electron states in agreement with Eqs. (I23l) - fl26l) ): 


hm G(X^,Xi;p+|7l) = 

T—>-00 

|£ Cj ^|—>-00 


1 /2hp+ 2^ _ ^ 

cr=l 




P+ > 0, 


hm G{Xa.X'yp4A) = 

i_ 


T'^-oo 

>-CX) 


hm G{Xa,X'y-\Pa\\A) 

T'^oo 

|aj^ 1^00 


hm G{Xa,X'y-\p+\\A) 

T^—oo 
|a:j^ |—>00 


(91) 

1 piiP±-m'^)\T'\/2hp+ ,2 _ 

— P+ > 0’ 


Awn? \T'\ 


cr=l 


1 /2h\p+\ ^2 ^_ 

— Y'/^ep°;Ux\A) 


Awh^ T' 


U-a- 


q,-a 


(92) 


V+ < 0, 


<7=1 


(93) 


1 e 


i{q'_^-'m?)\T\/2h\p+\ J? 




47rh2 |T| 


<7=1 


(94) 


where 'ip{x) = for a generic bispinor 'ijj{x) and where, consistently with the results in 

Section m the hrst-order terms in l/p+ in the pre-exponent proportional to the unit matrix 
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have been neglected (the four-momenta p'’^, q'^, and q’’^ are given in Eqs. (l8^ - fl8^ b 
Obviously, the in-and out-states in Eqs. fl9l]) - fl9^ differ from those in Eqs. fl2^ - fl26|) only 
by the 0-dependent phase factors exp(—ip+0) and exp( 2 p+ 0 ) for the positive- and negative- 
energy states, respectively, i.e., = exp(=Fip+0)-0^p('^(f^(X|y4). Also, it should 

be noticed that in Eqs. the quantity e = \/m? + is always a positive quantity 

and the appearance of the quantity \/^ instead that of 2|e| as in 8^ is only due to a 
difference in the normalization of the constant bispinors. Finally, it is also worth observing 
that, since in deriving the propagators we have neglected all quantum corrections in the 
exponents, the limiting procedure cannot reproduce the correcting terms proportional to 
the unit matrix in the pre-exponent in Eq. fll7p . as those arising, e.g., from the Van Vleck 
determinant. However, as we have mentioned in Section [TTl such terms can be neglected in 
the calculation of the probabilities of strong-held QED processes. 


IV. CONCLUSIONS AND OUTLOOK 


In conclusion, in the present paper we have continued and extended the investigation 


started in 


72| . and provided analytical tools for stuping strong-held QED processes 


m 


the presence of tightly focused laser beams. As in 72|, we have focused on the regime 
of interaction where the asymptotic energy of the electron is the largest dynamical energy 
scale in the problem. First, we have presented an equivalent but simpler expression of the 


wave functions found in 


72| . We have also generalized the expression of the wave functions. 


including the possibility that the electron also has a non-vanishing asymptotic transverse 
momentum. This is particularly important for the study of hrst-order strong-held QED 
processes, where the two leptons (e.g., initial and hnal electron or hnal electron and positron) 
do not propagate exactly along the same direction. Moreover, we have shown that the wave 
functions found here have essentially the same structure of the Volkov states and that the 
former can be practically obtained from the latter by means of a simple substitution of the 
transverse vector potential of the laser held. 

Moreover, we have now constructed the corresponding scalar and spinor propagators 
under the same approximations and we have verihed by means of an appropriate limiting 
procedure that the determined propagator is compatible with the corresponding wave func¬ 
tions. We have seen that the structure of the propagator is also essentially the same as the 
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Volkov propagator, with an effective transverse vector potential also containing derivatives 
of the focused laser field with respect to the transverse variables. This has allowed for the 
introduction of a local laser-induced mass correction also in the case of a focused laser field. 

In an upcoming publication we will apply the wave functions found here and in to 
investigate the first-order strong-held QED processes (nonlinear Compton scattering and 
nonlinear Breit-Wheeler pair production) including the tight focusing of the background 
laser held. 
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